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Abstract 

We describe the recent developments to extend the multi-parton in- 
teraction model of underlying events in Herwig-i-i- into the soft, non- 
perturbative, regime. This allows the program to describe also mini- 
mum bias collisions in which there is no hard interaction, for the first 
time. It is publicly available from versions 2.3 onwards and describes 
the Tevatron underlying event and minimum bias data. The extrapo- 
lations to the LHC nevertheless suffer considerable ambiguity, as we 
discuss. 

1 Introduction 

In this talk, we will summarize the development of a new model for the underlying event in Her- 
wig-i-i-, extending the previous perturbative multi-parton interaction (MPI) model down into the 
soft non-perturbative region. This allows minimum bias collisions to be simulated by Herwig-i-i- 
for the first time. 

We begin, though, by mentioning a few of the features that accompanied it in the re- 
lease of Herwig-i-i- [1] version 2.3 [2] in December 2008, which include NLO corrections in the 
POWHEG scheme for single W and Z production [3], and Higgs production [4]. Lepton-hadron 
scattering processes have been included for the first time. The simulation of physics beyond the 
standard model (BSM) has been extended to include a much wider range of 3-body decays and 
off-shell effects [5]. The treatment of baryon decays has been extended to match the sophistica- 
tion of meson and tau decays, including off-shell and form factor effects and spin correlations. 
Finally, in addition to the soft interactions discussed here, the MPI model has been extended 
to include the possibility of selecting additional scatters of arbitrary type, which can be impor- 
tant backgrounds to BSM signatures for which the single-scattering backgrounds are small, for 
example two like-sign Drell-Yan W productions [6]. 

The semi-hard MPI model was implemented in Herwig-i-i- version 2.1 [7]. It allows for 
the simulation of underlying events with perturbative scatters with pt > p™™ according to the 
standard QCD matrix elements with standard PDFs, dressed by parton showers that, in the ini- 
tial state, account for the modifications of the proton structure due to momentum and flavour 
conservation. It essentially re-implemented the existing Jimmy algorithm [8] that worked with 
the fortran HERWIG generator [9], but gave a significantly better description of the CDF data 
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on the underlying event [10], in part due to a more detailed global tuning [11]. However it was 
only able to describe the jet production part of the data, above about 20 GeV, and not the min- 
imum bias part, owing to a lack of soft scatters below p™™. A possible extension into the soft 
regime was first discussed in Ref. [12], but we have provided the first robust implementation 
of it, described in detail in Ref. [6]. It is somewhat complementary to the approach used in 
Pythia [13, 14], where the perturbative scatters are extended into the soft region through the use 
of a smooth non-perturbative modification. However, we make a stronger connection with infor- 
mation on total and elastic scattering cross sections, available through the eikonal formalism, to 
place constraints on our non-perturbative parameters [15]. 

In the remainder of this introduction. 
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we recap the basics of the eikonal model 
and recall the results of the perturbative MPI 
model that we had previously implemented in 
Herwig-i-i-, before showing how to extend it 
into the soft region. In Sect. |2] we discuss the 
constraints that can be placed on the model by 
the connection with hadronic scattering, and 
in Sect. [3] we show the predictions for final 
state properties. 

The starting point for the MPI model is 
the observation that the inclusive cross section 
for perturbative parton scattering may exceed 
the total hadron-hadron cross section. We 
show an example in Fig.[Tl with two of the to- 
tal cross section parameterizations we will be 
using. The origin of the steep rise in the par- 
tonic cross section is the proliferation of par- 
tons expected at small x. The excess of the 
partonic scattering cross section over the total 

cross section simply implies that there is on average more than one parton scattering per inelas- 
tic hadronic collision, n = crjet/cinei- Since the majority of scatters come from very small x 
partons, they consume relatively little energy and it is a good approximation to treat them as 
quasi-independent. 

From the optical theorem, one derives a relationship between the Fourier transform of the 
elastic amplitude a(b, s) and the inelastic cross section via the eikonal function, x(b, s), 
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Fig. 1: Total cross sections (black) in the two parameteri- 
zations of Donnachie and Landslioff [16, 17]. In blue the 
QCD jet production cross section above 2 GeV is shown. 
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One can construct a QCD prediction for the eikonal function by assuming that multiple scatters 
are independent, and that the partons that participate in them are distributed across the face of the 
hadron with some impact parameter distribution G{h) that is independent of their longitudinal 
momentum. 



XQCD(b, s) = i ^(b) aZ,^is), A{h) = J dV G(b') G(b - b') 



(2) 



where cr\f^^^ is the inclusive partonic scattering cross section, which is given by the conventional 
perturbative calculation. 

In the original Jimmy model and its Herwig++ reimplementation, these formulae are im- 
plemented in a straightforward way, with the hard cross section defined by a strict cut, pt > p™™ 
and the matter distribution given by the Fourier transform of the electromagnetic form factor. 



with, to reflect the fact that the distribution of soft partons might not be the same as that of 
electromagnetic charge, /u^ considered to be a free parameter and not fixed to its electromagnetic 
value 0.71 GeV^. Compared to a Gaussian of the same width, this distribution has both a stronger 
peak and a broader tail so it is somewhat similar to the double-Gaussian form used in Pythia [18]. 
In Ref. [15], we explicitly showed that the two result in similar distributions, if their widths are 
fixed to be equal, except very far out in the tails, /i^ and pf™ are the main adjustable parameters 
of the model and, allowing them to vary freely, one can get a good description of the CDF 
underlying event data, as shown in Fig.[2l The choice of parton distribution function can also be 
seen to have a small but significant effect. 

The main shortcoming of this model is that it does not contain soft scatters and hence 
cannot describe very low pt jet production or minimum bias collisions. In Ref. [12] it was 
proposed to remedy this, by extending the concept of independent partonic scatters right down 
into the infrared region. One can therefore write the eikonal function as the incoherent sum of 
the QCD component we already computed and a soft component. 



Xtot(b,s) = XQCD(b,s) + Xsoft(b,s) = ^(^A{h)aZrd{s)+A,oit{h)aZ{t{s)), (4) 



where (T™£^ is an unknown partonic soft scattering cross section. As a first simplest model, we 
assume that the matter distributions are the same, Asoft(b) = A{h), although we relax this 
condition later. By taking the eikonal approach seriously, we can trade the unknown soft cross 
section for the unknown total hadronic cross section. 



Knowing the total cross section, for a given matter distribution and hard cross section (implied 
by p™™ and the PDF choice) the soft cross section is then determined. In order to make pre- 
dictions for energies higher than the Tevatron, we consider three predictions of the total cross 
section: 1) the standard Donnachie-Landshoff parameterization [16]; 2) the latter for the energy 
dependence but with the normalization fixed by the CDF measurement [21]; and 3) the newer 
Donnachie-Landshoff model with a hard component [17]. Of course once we have an experi- 
mental measurement from the LHC we would use that for our predictions. In this way, our simple 
hard-i-soft model has no more free parameters than our hard model and we can tune //^ and p™™. 
Before doing this, we present the results of Ref. [15], in which we considered the theoretical 
constraints that could be put on these parameters. 
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Fig. 2: Multiplicity and Pt"™ in the transverse region. CDF data are shown as black circles, Herwig++ without 
MPI as magenta dots, with MPI using MRST [19] PDFs as solid red and with CTEQ6L [20] as cyan dashed. The 
lower plot shows the statistical significance of the disagreement between the Monte Carlo predictions and the data. 
The legend on the upper plot shows the total for all observables, whereas the lower plot for each observable has its 
values. 

2 Analytical constraints 
2.1 Simple model 

Within our model we want (T™f^^ to correspond to a physical cross section. It must therefore be 
positive. This therefore places constraints on the /i^-p™™ plane: a lower bound on p™™ for a 
given value of jj?. These are shown for the Tevatron on the left-hand side of Fig. [3] as the sohd 




Fig. 3: Left: The parameter space of the simple eikonal model at the Tevatron. The solid curves come from u"a{t > 
for three different PDF sets. The horizontal lines come from&ci = 16.98 ± 0.25 GeV~^ [21,23]. The excluded 
regions are shaded. The dashed lines indicate the preferred parameter ranges from the fit to Tevatron final-state 
data [11]. Right: The equivalent plot for the LHC. The additional (dashed) constraints come from requiring the total 
number of scatters to be less than 10. 



lines for three different PDF sets: the two shown previously and IVIRST LO* [22]. Since in the 
eikonal model the total and inelastic cross sections are related to the elastic one, we can also 
place constraints from the elastic slope parameter, which has been measured by CDF [21,23]: 
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(17 ± 0.25) GeV" 



(6) 



This rather precise measurement directly constrains fi^ in our simple model and rules out all but 
a very narrow strip of the parameter space. Finally, we consider the parameter space of the fit 
to final-state data. Although there is a preferred point in the parameter space, the tuning of both 
the hard-only model [11] and the hard-i-soft model shown below indicates a strong correlation 
between the two parameters and there is a broad region of acceptable parameter values, which 
we show in Fig. [3]by the region edged by red bands. Between the different constraints we have 
only a very small allowed region of parameter space. 

At the LHC the picture is similar, although the constraint dg^f^ > is considerably more 
restrictive (note the difference in range of the x axes of the two plots). Different models predict b^i 
in the range 19 to 22 GeV^^ translating into a slightly wider horizontal band. Finally, although 
we do not have final-state data to compare to, in order to simulate self-consistent final states at 
all we find that we must prevent the multiplicity of scatters becoming too high. While precisely 
where we place this cut is arbitrary, we indicate it by shading the region in which the mean 
number of scatters is greater than 10. This plot is shown for the central of the three total LHC 
cross section predictions we consider - it is qualitatively similar for the other two, although the 
different constraints move somewhat. 



Comparing the two plots in Fig. [3l we come to the reahzation that, from these theoretical 
constraints together with the fit to the Tevatron data, we can already rule out the possibility that 
the parameters of this simple model are energy-independent - there is no region of the plot that 
is allowed at both energies. 

While it could be that the parameters of the MPI model are in fact energy dependent, as 
advocated by the PYTHIA authors [24], we prefer to let the LHC data decide, by proposing a 
model that is flexible enough to allow energy-independent or -dependent parameters. The sim- 
plest generalization of the above model that achieves this is actually well physically motivated, 
and we call it the hot-spot model. 

2.2 Hot-Spot model 

The simple model has other shortcomings, beyond our aesthetic preference to allow the possi- 
bility of energy-independent parameters. The values of 0"™^^^ extracted from the predictions of 
o"tot [15], have rather strange energy dependence, being quite sensitive to precise details of the 
matter distribution, parameter choice, cross section prediction and PDF set and, in most cases, 
having a steeply rising dependence on energy, much steeper than one would like to imagine 
for a purely soft cross section. Moreover, the value of fi^ extracted from b^i is in contradic- 
tion with that extracted from CDF's measurement [25] of double-parton scattering, which yields 
= 3.0 ± 0.5 GeV^. 

All of these shortcomings can be circumvented by allowing the matter distribution to be 
different for soft and hard scatters. As a next simplest model, we keep the same form for each, 
but allow the /i^ values to be different. We again fix the additional free parameter, this time to a 
fixed value of 6ci- That is, once atot and 6ei are measured at some energy, the non-perturbative 
parameters of our model, (T™£j and fi^^^^ are known. Since it will turn out that our preferred value 
of is significantly larger than the extracted value of ^g^f^, we call this a hot-spot model: soft 
partons have a relatively broad distribution, actually similar to the electromagnetic form factor, 
while semi-hard partons (typically still small x, but probed at momentum scales above ™) are 
concentrated into smaller denser regions within the proton. 

Having used one constraint to fix an additional parameter, there is only one constraint 
in the parameter space, shown in Fig. |4] for the Tevatron and LHC. The model has much more 
freedom than the simple one, with much of the parameter space allowed, and with ample overlap 
between the allowed regions at the two energies. 

Another nice feature of this model is the energy-dependence of cr^j^ it implies, shown in 
Fig. m At least for the standard Donnachie-Landshoff energy dependence, it corresponds to a 
very slow increase, almost constant, in-keeping with one's expectations of a soft cross section. 

3 Final states 

We have implemented this model into Herwig-i-i-. There are many additional details that we 
do not go into here [6], but wherever possible, the treatment of soft scatters is kept as similar 
as possible to that of semi-hard scatters, to make for a smooth matching. In particular, for the 
transverse momentum dependence, we make the distribution of a Gaussian centred on zero, 
whose integral over the range zero to pf™ is given by a"^j^ and whose width is adjusted such 




Fig. 4: Parameter space of the improved eikonal model for the Tevatron (left) and LHC (right). The solid curves 
impose a minimum allowed value of /i^, for a given value of p™'" by requiring a valid description of crtot and bc\ with 
positive (Tgoft. The excluded regions are shaded. We used the MRST 2001 LO [19] PDFs for these plots. 



that da/dpt is continuous at p™™ is therefore seen to be not a cutoff, as it is in the Jimmy 
model, but a matching scale, where the model makes a relatively smooth transition between 
perturbative and non-perturbative treatments of the same phenomena, in a similar spirit to the 
model of Ref. [26] for transverse momentum in initial-state radiation. 

The model actually exhibits a curious feature in its pt dependence, first observed in Ref. [12]. 
With the typical pai^ameter values that are preferred by the data, da/dpt is large enough, and Ug^f^ 
small enough, that the soft distribution is not actually a Gaussian but an inverted Gaussian: its 
width-squared parameter is negative. The result is that the transverse momentum of scatters is 
dominated by the region around p™™ and not by the truly non-perturbative region pt —>■ 0. This 
adds to the self-consistency of the model, justifying the use of an independent partonic scattering 
picture even for soft non-perturbative collisions. 

With the model in hand, we can repeat the tune to the CDF data on the underlying event. 
Unlike with the semi-hard model, we now fit the data right down to zero leading jet momentum. 
The result is shown in Fig. [6l which is qualitatively similar to the one for the semi-hard model. 
The description of the data in the transverse region is shown in Fig. |7J It can be seen to be 
reasonable in the lower transverse momentum region, although certainly still not as good as at 
higher transverse momenta. 

The discrepancy in the lowest few bins may be related to another deficiency of our model. 
According to the eikonal model, the inelastic cross section should include all final states that 
are not exactly elastic, while our simulation of them generates only non-diffractive events in 
which colour is exchanged between the two protons and hence a significant number of final- 
state hadrons are produced. While single-diffractive-dissociation events would not be triggered 
on experimentally, double-diffractive-dissociation events, in which both protons break up but do 
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Fig. 5: (Jgoft as a function of energy. Each of the three Fig. 6: Contour plots for the per degree of freedom 

different curves shows the soft cross section that would for the fit to the CDF underlying event data. The cross 
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cross section is used. Curves that do not reach out to 30 area consists of parameter choices where the elastic t- 
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reproduce atot and b^i correctly at these energies. simultaneously. 



not exchange colour across the central region of the event, would, and would lead to extremely 
quiet events with low leading jet pt and low central multiplicity, which are not present in our 
sample. In Ref. [6] we have checked that these bins are not pulling our tune significantly by 
repeating it without them. The overall chi-squared is significantly smaller, but the best fit point 
and chi-squared contours are similar. 

4 Conclusions 

We have reviewed the basis of the semi-hard MPI model that we previously implemented in 
Herwig-i-i-, and motivated its extension to a soft component. Through the connection with the 
total and elastic cross sections provided by the eikonal model and optical theorem, we have 
placed significant constraints on the simplest soft model. We have shown that these constraints 
can be relaxed by invoking a hot-spot model in which the spatial distributions of soft and semi- 
hard partons are different. Finally, we have implemented this model and shown that it gives a 
reasonable description of the minimum bias data, for the first time in Herwig-i-i-. Nevertheless, 
there is still room for improvement, particularly in the very low pt region and several avenues for 
further study present themselves, not least the diffractive component already mentioned, and the 
role of colour correlations, which were argued to be very important in Ref. [14], but which seem 
to be less so in the current Herwig-i-i- implementation [6]. 

Despite the successful description of Tevatron data, the extrapolation to the LHC suffers 
from considerable uncertainty. The unknown value of the total cross section, which determines 





Fig. 7: Multiplicity and p^"™ in the transverse region. CDF data are shown as black circles. The histograms show 
Herwig++ with the improved model for semi-hard and soft additional scatters using the MRST 2001 LO [19] PDFs 
for three different parameter sets. The lower plot shows the ratio Monte Carlo to data and the data error band. The 
legend shows the total for all observables. 

the non-perturbative parameters in our model, plays a crucial role, but even once this and the 
elastic slope parameter have been directly measured, the region of allowed parameter space is still 
large. Although we prefer a model in which the parameters are energy independent, ultimately 
only data will tell us whether this is the case. Finally, even once the underlying event data have 
been measured, the parameters will not be fully tied down, due to their entanglement with the 
PDFs. We eagerly await the LHC data to guide us. 
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